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The unusual electronic properties of single-layer graphene[lj make it a promis- 
ing material system for fundamental advances in physics, and an attractive plat- 
form for new device technologies. Graphene's spin transport properties are 
expected to be particularly interesting, with predictions for extremely long co- 
herence times and intrinsic spin-polarized states at zero field [2, El HI [5]. In 
order to test such predictions, it is necessary to measure the spin polarization 
of electrical currents in graphene. Here, we resolve spin transport directly from 
conductance features that are caused by quantum interference. These features 
split visibly in an in-plane magnetic field, similar to Zeeman splitting in atomic 
and quantum dot systemspj], \7\. The spin-polarized conductance features that 
are the subject of this work may, in the future, lead to the development of 
graphene devices incorporating interference-based spin filters. 

Spin currents have recently been measured in graphene using ferromagnetic contacts, an 
advance that shows promise for room-temperature graphene spintronics[Hl El HD]- Unfor- 
tunately, many of the other tools that have been developed for studying spin-dependent 
phenomena in solid-state nanostructures (optical excitation, spin-resolved Coulomb block- 
ade, etc.) have proven difficult to implement in graphene. The Zeeman-split conductance 
fluctuations that are the subject of this paper provide a way to detect spin-polarized currents 
with a straightforward electrical measurement. This is the first report, to our knowledge, of 
visible spin-splitting of conductance fluctuations in any material. 

Flakes of graphene were deposited using well-established mechanical exfoliation 
techniquesp] onto silicon wafers with a ~300nm surface oxide, then contacted with Cr/Au 
leads (Fig. la). The flake used for this experiment was selected optically and confirmed to 
be single layer by quantum Hall measurements (see supplement) showing the characteristic 
(n+1/2) Landau quantization[lJ. The two-terminal conductance measurements described 
below were performed in a dilution refrigerator with a base temperature of 20mK. The dis- 
tance between the contacts was 3/zm; the phase coherence length was estimated to equal or 
exceed sample dimensions. A two-axis magnet provided fields up to S||=12T in the plane 
of the graphene and fine control of the out-of-plane component up to Sx=120mT (Fig. la). 

Carrier density n was controlled capacitively with a voltage Vq on the backgate, inducing 
a density n = aVc in the graphene sheet. The proportionality constant, a = 8.05±0.05x 10 10 
cm~ 2 /V, was extracted from the quantum Hall data. A Dirac resistance peak was observed 
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when Vg tuned the carrier density to the charge neutrality point at Vo~lV. The conduc- 
tance, G, was calculated after taking into account an estimated contact resistance of 3.2kQ. 
Conductance fluctuations, SG, were visible superimposed on the Dirac lineshape (Fig. lb), 
and were extracted by subtracting an ensemble-averaged background (see supplement). 

Graphene, like any metallic or semiconducting system, exhibits 'universal' conductance 
fluctuations when transport is phase coherent [TT| [T2l IT3] . These aperiodic fluctuations 
arise from quantum interference between all possible paths that charge carriers take as they 
traverse the device. The interference pattern is randomized by small changes in the relative 
phase between trajectories, which depends on the Aharonov-Bohm magnetic flux through 
the sample as well as on the Fermi wavevector|llj. The conductance includes contributions 
from spin-up and spin-down carriers: G = + G±. In the absence of spin-orbit and 
many-body interactions, G-m depends on the spin only through the Fermi wavevector, and 
therefore the density, of the respective carrier population. 

Spin-up and spin-down carriers have the same density at zero field, n^=n^. In a simple 
picture, then, zero-field conductance fluctuations include identical contributions from both 
populations (Fig. lc). Applying a magnetic field partially polarizes the graphene - it induces 
a difference between spin- up and spin-down carrier densities. The total density n = n^+ni = 
aVc is set by the backgate voltage; the difference in the densities n^—ni = ^^gusB is set 
by the magnetic field, where 4^ is the density of states and g^B is the Zeeman energy. 
Together these yield 

» - Q (v + V ^ - 1 dn nn R (1) 

n m-^[Vc±—), v G -- Te g^B (1) 

As a result, spin-up and spin-down conductance contributions at finite field are offset in gate 
voltage by V^ ffset (Fig. Id), leading to Zeeman splitting of interference features in a gate 
voltage trace G{Vq) (Fig. le). 

Zeeman splitting of graphene conductance fluctuations can be seen in Figs. IF and 1G. 
The 'V shapes in the data correspond to spin-resolved conductance: left-moving (right- 
moving) arms reflect interference for a particular density of spin-down (spin-up) carriers. 
Although conductance fluctuations have been observed for several decades in a wide variety 
of materials [HI [T3], a direct measurement of their Zeeman splitting has never before been 
reported. We speculate that this effect is visible in graphene due in part to its atomic-scale 
flatness and low density of states. 
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A statistical analysis of the spin-split offset at particular values of B» was obtained using 
an ensemble of traces, 8G(Vq), collected for a range of out-of-plane fields, 3mT<i?x<120mT 
(•C-B||). The offset is barely visible in the raw conductance data (Fig. 2a), but can be seen 
more clearly by computing an autocorrelation C[sg\(AVq) = (5G(Vg)SG(Vg + AVq))- Side- 
peaks in the autocorrelation at AVg = ±VQ ffset (Fig. 2b) reflect similar features offset in the 
conductance traces (Fig. 2a) due to the difference in spin-up and spin-down carrier densities. 

The side-peaks shifted outward with increasing Bn (Fig. 3a), while the height of the 
central peak (C[sc\(0), the variance of the fluctuations) dropped by a factor of two (Fig. 3b). 
The drop in C\sq\(0) reflects a suppression of conductance fluctuations when the Zeeman 
energy, g^B, exceeds energy broadening due to temperature or dephasing. This effect has 
been used to characterize spin degeneracy in other systems [T4" l [T5l fT6] . where splitting could 
not be observed directly. The suppression factor of ~2 indicates that spin degeneracy was 
intact at zero field, after taking into account broken time reversal symmetry due to the small 
out-of-plane field. The shift of the side-peak location, V^ ffset , with B» was expected from 
Eq. 1 — a statistical confirmation of the linear spin-splitting seen in Figs. IF and 1G. 

The density of states can be extracted directly from VQ ffset at a given field using g=2 
for graphene (see supplement) and the value for a determined from quantum Hall mea- 
surements (Eq. 1). VQ ffset (Vcj) was recorded over the full gate voltage range by computing 
the autocorrelation C[sc\(AVq) within a sliding window in Vq (Fig. 3c). The resulting 
square-root lineshape can be compared to 4^ expected from graphene's dispersion relation: 
^(Vg) = 2 ^{^ \J\Vg — Vo\, with a Fermi velocity vp that would be independent of density 
for an ideal Dirac band structure. 

A quantitative analysis of VQ ffset (VG) gives fp>1.0xl0 6 m/s over the range of densities 
reported here (Fig. 3d), consistent with measurements using a variety other techniques [TTl 
[HI US], but considerably larger than the expected Vp ~ 0.8xl0 6 m/s from non-interacting 
band structure calculations [4J. Enhancements in v-p have been predicted due to many-body 
effects [201 [2Tj . which become particularly important at low densities. Indeed, larger values 
of v-f are observed in Fig. 3d at lower densities, similar to a report of infrared spectroscopy 
measurements [18]. An additional deviation from the ideal Dirac density of states is seen at 
charge neutrality (from Vq ~ 0-2V): rather than dropping to zero, a finite value of VQ ffset 
is visible (Fig. 3c, supplement). This minimum density of states, 4 x 10 12 cm" 2 /eV, is 
similar to that observed in scanning single-electron transistor measurements [19J, attributed 
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to broadening from disorder. 

In order to observe V's like those shown in Figs. IF and 1G, interference features must 
shift, but not otherwise change, as a function of in-plane field. In other words, the primary 
influence of the magnetic field on the interference must be through its effect on the densities 
of spin-up and spin-down carriers, rather than a change in Aharonov-Bohm (AB) flux. 
Fields of a Tesla or more were required to differentiate the two carrier densities, but the AB 
flux introduced by a few milliTesla in B± was sufficient to alter the observed interference 
pattern. The drastically different field scales for spin and AB effects demanded precise 
alignment (within 0.05°) of the B» axis to the plane of the graphene flake. The alignment of 
B\\ was monitored using weak localization (WL), the coherent enhancement of backscattering 
associated with time-reveral symmetry at B i =0[TT] . Slight mis- alignment led to a shift in 
the location of the WL conductance dip as B» was raised, and was corrected by offsetting 
B ± . 

In addition to shifting, the WL dip also decreased in magnitude with Bn, by a factor of 
two at IT and below detectable levels above 4T (Fig. 4ab). The complete collapse of the WL 
dip, over a field range where the variance of conductance fluctuations decreases only by a 
factor of two, implies that time-reversal symmetry is broken even by an in-plane field [22| 123]. 
The disappearance of symmetry in g(B±) <-> g(—B±) at finite B\\ (Fig. 4cd) provided further 
evidence of time-reversal symmetry breaking|22j. Similar phenomena have been observed 
in semiconductor 2D electron gases, and have been associated with finite thickness and 
nanometer-scale undulations in the 2DEG[23J, allowing B\\ to thread an AB flux through 
the conductor. The -By-scale of WL collapse in graphene corresponds to an effective thickness 
of ~lnm, in agreement with previous measurements of the intrinsic ripple size in graphene 
sheets [231 125]. 

We acknowledge V. F'alko, D. Goldhaber-Gordon, H. Heersche, T. Ihn, G. Kamps, C- 
Y. Lo, A. Morpurgo, and K. Todd for helpful discussions. Graphite provided by S. Fain and 
D. Cobden. MBL acknowledges a PGS-D from NSERC; work funded by CIFAR, CFI, and 
NSERC. 
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Fig. 1. Zeeman effect in graphene. (a) Schematic of single- layer graphene device, show- 
ing orientation of magnetic fields B\\ and B±. (b) Two-terminal conductance, G, measured 
as a function of backgate voltage Vq at 20mK, showing the Dirac conductance minimum 
(i?^=5|l=20mT). Inset: a sample of conductance fluctuations over a narrow range of gate 
voltage. (c,d,e) Diagrams of spin-splitting effect: a simulated conductance trace at zero 
field (c); two offset contributions from spin- up and spin-down at finite field (d); contributions 
add together in the total conductance (e), giving spin-split features. (f,g) Two examples of 
conductance fluctuations spin-splitting in an in-plane field, taken over two ranges of gate 
voltage (-B_i_=80mT). For both images, a smooth background was subtracted to highlight 
fluctuations. The offset (indicated by dashed lines as guides to the eye) due to spin is larger 
around -6V compared with -2V because density of states increases with density: the dashed 
lines correspond to densities of states of (f) 12xl0 12 and (g) 7x 10 12 cm~ 2 /eV. 
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Fig. 2. Autocorrelation extraction of spin splitting, (a) Features are doubled in an 
image of conductance fluctuations at B»—8T. (b) Autocorrelations of the windows indicated 
in (A) show that the spin-split offset (the location of the autocorrelation side-peak) changes 
from V^ ffset =0.22V at V G =-16V, to V^ ffsct =0.19V at Vb=-12V. The autocorrelation function 
is symmetric, so side-peaks appear at ±VQ ffsct . 
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Fig. 3. Magnetic field- and gate-dependence of spin splitting, (a) Autocorrela- 
tions, C[sg](AVg), of fluctuations in the range Vg=-5...-4.4V, averaged over an ensemble 
of traces for 3mT<5x<120mT, demonstrate the in-plane field dependence of side-peak 
location, V^ ffset . The dashed red line is fit to VQ Sset , and corresponds to a density of 
states 1.04xl0 13 cm~ 2 /eV. (b) Variance C[sg\(0) (filled squares) and autocorrelation side- 
peak height C[SG]{VQ Gsct ) (open squares) from (a). Variance drops to the expected 50% of 
its zero-field value, but the side-peaks fall below the expected 25%. The out-of-plane field 
-Bj_>3mT was large enough to break time reversal symmetry even at JBii = 0. (c) Auto- 
correlations computed for a sliding 0.75V- wide window in Vq (c.f. Fig. 2), and averaged 
over an ensemble of traces for 3mT<i?x<120mT, show how the density of states depends 
on density, (d) The Fermi velocity vp extracted from side-peak locations, VQ Sset , in (c) is 
enhanced near the Dirac point. Error bars are primarily due to uncertainty in the charge 
neutrality point Vb=1.0±0.3V. 
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Fig. 4. Loss of weak localization and symmetry due to in-plane field, (a) Aver- 
aging out the conductance fluctuations over a range Vg=-5...-3V, a weak localization dip 
in conductance is present at B=0 but disappears for increasing B\\ or B±. This indicates 
a breaking of time-reversal symmetry from either component, but at very different field- 
scales. Curves at different B\\ are vertically offset to align G(B j_>10mT) (see supplement), 
(b) Negative magnetoconductance for 5x>10mT is unaffected by B». Box encloses range 
from (a). (c,d) Conductance is symmetric in B± at B\\=0 (c) but not at Bn—2T (d). These 
data provide further evidence of time-reveral symmetry breaking by the in-plane field, ex- 
cluding (for example) that the disappearance of weak localization in (a) and (b) is due to 
decoherence. 
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Supplementary information 

[Device Fabrication] Graphene was deposited on SiO^/Si wafers made from highly- 
arsenic-doped silicon with resistivity < 0.007 f2-cm, thermally oxidized by the supplier to 
a uniform depth of 300nm±15%. Before graphene deposition, the oxide was thinned using 
a CF4/O2 plasma to change the colour from the initial teal to purple, ~ 265nm, to aid in 
the identification of deposited graphene. After thinning, wirebonding pads and alignment 
features were laid down by standard optical lithography techniques. The surface was cleaned 
with a standard SC-l/SC-2 etch immediately before graphene deposition. 

Natural graphite flakes were cleaved repeatedly by Scotch tape, then pressed lightly on 
the chip surface. After 10 minutes, the tape was slowly peeled off the chip, then the chip 
was placed in hot acetone for 2 minutes to remove nearly all the tape residue, followed 
by isopropanol rinse and drying. Graphene flakes were identified by their colour in an 
optical microscope, where colour shifts were observed in roughly equal steps for 0,1,2,3 
layers. After identification, flakes were contacted by thermally evaporated 5nm Cr / 50nm 
Au metallization, patterned using an electron-beam lithography process at 20keV with a 
200K/950K PMMA bilayer resist (Supp. Fig. [TJ. 

Several single- layer flakes were evaluated based on the width and offset [26] of the Dirac 
resistance peak R(Vg). Typical widths were 5-10V, centered around a voltage often less 
than 5V. These characteristics were achieved only with the SC-l/SC-2 etch step mentioned 
above — without this step, flakes were ~50V wide and centered at >30V. The flake measured 
in this paper had a narrow but slightly asymmetric Dirac peak, with local maxima at +0.5V 
and +1.5V (Supp. Fig. [2}. 

[4K and 20mK cooldowns] The flake was first measured at 4K in a quantum Hall 
configuration. The characteristic n + 1/2 quantization was observed, confirming that the 
flake was single layer graphene pQ (Supp. Fig. [3]). The gate dependence of Shubnikov-de-Haas 
oscillations allowed an accurate determination of the induced carrier density coefficient a 
(Supp. Fig. [3]). The measured value of a was crucial for quantifying the density of states 
and Fermi velocity throughout the paper. The linear relationship n = oVq is a result of the 
geometrical capacitance c g = ae = (esi02/thickness) between the backgate and graphene. 
In general, with finite densities of states an additional nonlinear 'quantum capacitance' 
c q = e 2 ^ is present in series with the geometrical capacitance [27J. Since the gate oxide 
was thick, the correction to total capacitance was only 2% for the lowest density of states 
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observed. Hence, variations of a were ignored in our analysis. 

During the transfer from the 4K probe to the dilution refrigerator one of the contacts to 
the graphene was broken, preventing a 4- wire resistance measurement. As a result, the data 
presented in the main paper are two-terminal measurements. 

[Conductance Measurement] Dilution refrigerator measurements were performed us- 
ing lOnA lock-in excitation at 77.92 Hz. A contact resistance of 3200 Q was estimated by 
extrapolating the Dirac peak resistance to large negative gate voltages (Supp. Fig. [2]); this 
resistance was subtracted before calculating the conductance. The Dirac peak extrapolated 
to 4200 Q for positive gate voltages, perhaps due to a diode effect in the contacts. This Vg- 
dependent contact resistance does not affect our analysis, as it is insensitive to the absolute 
amplitude of SG fluctuations. Mobility was estimated using the conductance data in Fig. lb 
and assuming a two-square geometry, giving 4000 cm 2 V _1 s _1 for high carrier densities. 

[Field orientation] The two-axis magnet used for dilution refrigerator measurements 
had a primary solenoid aligned to within ~ 0.4° of the plane of the sample, controllable up 
to £>||=12T, and a split-coil magnet providing a field up to B± = 120mT perpendicular to 
the sample plane. The separate effect of the two magnetic field components can be seen in 
Supp. Fig. [4J showing symmetry in total field reversal B <-> —B as required by reciprocity 
relations |28l 129] for a two-terminal measurement. 

The B± magnet was used to correct for the 0.4° error in B» alignment, measured precisely 
by the shift in the weak localization conductance (WL) dip and symmetry center. Although 
WL could not be measured accurately for B» > 3T, the B± correction above that field 
assumed the same mis-alignment angle as measured at lower field. Throughout the paper, 
the stated B± is the compensated field, as experienced by the graphene. Although the V's 
from Figs. If and lg were sensitive to small errors in B\\ alignment, the autocorrelations 
of gate voltage traces were insensitive because they were computed for a specific Bn, and 
averaged over a range of fields 3mT< B± < 120 mT. 

[Background subtraction] In order to extract the fluctuations SG = G — (G) from the 
conductance data G, it was necessary first to evaluate (G), the conductance background. 
Background subtraction was an important step in the calculation of autocorrelation func- 
tions, and had to be done carefully to avoid introducing artifacts. Two different background 
extraction protocols were used, for data sets covering very different parameter ranges, and 
are described below. In all cases, the details of background subtraction influenced the SG 
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variance and clarity of the autocorrelations but did not affect the position of the sidepeaks, 
that is, the background did not affect the extracted V~Q ffset . 

The data in Figs. 2 and 3c were at B\\ = 8T, but covered a wide range of Vq- 
(Supp. Fig. [5]). For these data sets, a linear function (G)(B±) was fit at each value of 
Vg and subtracted off. Since the range of B± encompasses many conductance fluctuations, 
this procedure does not remove the primary signal, but it does remove larger gate- and 
field-scale fluctuations that would otherwise overwhelm the autocorrelation. 

For Fig. 3a, which covered only a narrow range of densities but a wide range of in- 
plane fields, the average effect on conductance of each of the three control parameters — 
(G)(Vg), (G)(B±), and (G)(B\\) — was calculated by averaging fluctuations over the other 
two parameters. For example, (G)(B\\) was calculated by averaging over both Vg and B±, 
denoted (G)v g ,b ± (B\\). Then the three-dimensional function (G)(Vq, B±, B») was calculated 
by adding the three effects together. The subtracted background was (G)(Vq, B±, B\\) = 



[Graphene g-factor] Graphene is expected to have a g-factor close to that of the free 
electron, but to our knowledge graphene's g-factor has been measured only in the quantum 
Hall regime|30j, where Landau quantization can lead to significant modifications [31]. On 
the other hand, extensive spin resonance measurements in graphite have found =2.0029 
for in-plane fields [32J [33J, [33] - We therefore approximate g=2 throughout the paper. 

[Density of states and Fermi velocity] The density of states for any 2D isotropic 
dispersion relation e{k) = e (y/k% + k^j with degeneracy D can be written as ^ = ^ /^| = 
y^Dn/n/ 1|. Defining Fermi velocity as Vp = |4§, this gives 



with graphene's D = 2 x 2 from valley and spin degeneracy. The Fermi velocities reported 
in Fig. 3d were calculated using: 

• the density of states ^ = - B VQ ffset extracted from the spin- splitting offset Vg ffset at 



• the density n = a(Va — Vq) determined by the gate voltage Vq away from charge 
neutrality Vq. 



(G)b ± ,bAV g ) + (G) b «,Vg(B±) + (G)v g ,b ± ( b \\) ~ 2 (G)v Q ,B ± ,B l 




B\\ = 8T (Fig. 3c, Supp. Fig. pj), and 
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Uncertainty in Vq was the primary source of the uncertainty in v-p in Fig. 3d. As seen in 
Supp. Fig. [6j a more positive value for V implies a reduced Vp for electrons (enhanced for 
holes); a more negative value does the opposite. The value used in the paper, V = 1.0±0.3V, 
was chosen to reflect the position, width, and asymmetry of the Dirac peak (Supp. Fig. [2]). 
Beyond this range, the extracted values for vf were highly asymmetric for electrons and 
holes at low density. 

[Deviations from ideal spin splitting] Some experimental observations point to a 
more complicated spin-splitting behavior than simple picture of two offset densities described 
in the text, indicating that interference for spin- up and spin-down carriers is somewhat 
different. First, the V's from Fig. If and lg in the main text contain more features than 
would be expected simply from offset densities. For comparison, the data from Fig. lg is 
replotted next to a simulation of "ideal" splitting based on the B» = trace from Fig. lg 
(Supp. Fig. [7]). Some features move either left or right in gate voltage with B\\ but do 
not split, that is, some left or right moving features do not appear to have an oppositely- 
moving counterpart. One interpretation of this observation would be that the conductance 
fluctuations were partially spin-resolved even at zero field, but further work is required to 
explore this possibility. 

Second, in Supp. Fig. [6] the side-peak amplitude in the autocorrelations is consistently 
smaller compared to the central peak than the 50% that would be expected for two offset 
copies of random fluctuations (as in Supp. Fig. [7j). One explanation would be that a shift in 
gate voltage changes not only the density but also the details of the scattering. One might 
expect from this that the relative height of the side-peak would fall off with field, as V^ set 
gets larger and larger and the offset features appear farther and farther from each other 
in gate voltage. In fact, the opposite behavior is seen in the data: the side-peak is higher 
at very high fields in Fig. 3b. As mentioned above, slight mis-alignment in B\\ would not 
explain the too-small side-peaks, which are taken from autocorrelations at a fixed value of 
B\\- 

[Finite density of states at charge neutrality] The spin-splitting offset does not 
drop to zero at the charge neutrality point as expected from a linear dispersion relation. 
Instead, a closer examination of the Vq = — 1...3V region from Fig. 3c (Supp. Fig. [8]) 
reveals a density of states 4^ = 3.5 ± 0.4 x 10 9 cm _2 /meV that is roughly independent of 
gate voltage. This consistent density of states is seen in both the correlation function and 
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the 'V's in the raw data. In scanning single-electron transistor measurements [19], the inverse 
compressibility (reciprocal of the density of states) was observed to reach a maximum value 
of about 3 xlO -10 cm 2 meV at the Dirac point, in agreement with our measurement. 

[Ripples] An order-of-magnitude estimate of the ripple size can be made by comparing 
the in-plane and out-of-plane magnetic field scales required to suppress weak localization. 
As seen in Fig. 4, the weak localization dip was supressed by 50% for B± of about 0.5 mT, 
corresponding to a a/ (/i/e)/0.5mT = 3/zm lateral phase-coherence length. In contrast, an 
in-plane field of approximately B\\ ~ 1.5T was required to reach 50% suppression, indicating 
that the effective vertical spread of backscattering trajectories (which have 3/im lateral 
length) is (^i/e)/(1.5T x 3/zm) = lnm. 

[Negative magnetoconductance from in-plane field] A significant decrease in the 
conductance for larger B\\ was observed over the full range of gate voltage, for any B± 
(Supp. Fig. |9j). The fractional change was -10% to -20% for B\\ up to 10T. An in-plane 
magnetoconductance has been predicted due to reduced screening in a spin-polarized system 
[35] but the observed effect was three orders of magnitude larger than predicted. The origins 
of this in-plane magnetoconductance are not yet understood. 
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Supplementary Figure 1: Optical microscope images of the graphene flake on SiC>2 before and 
after microfabrication. Scale bars: 3/im. Arrows indicate the two contacts used in the experiment. 
The other wires did not make contact at low temperature. 
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Supplementary Figure 2: Dirac peak in resistance from 20mK cooldown. The Dirac peak 
extrapolates to a different contact resistance for highly negative (blue, ~3200r2) or highly positive 
gate voltage (red, ~4200S7). Inset: A slight splitting in the peak is observed, possibly due to surface 
contamination of the graphene. 
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Supplementary Figure 3: Quantum Hall effect, (a) Longitudinal resistance as a function of 
gate and perpendicular magnetic field. Colour scale is chosen to emphasize Rxx minima due to SdH 
oscillations. The fan of lines corresponds to a=8.05xl0 10 cm~ 2 /V and charge neutrality at +0.4 
V. (b, c) Transverse hall resistance, showing the characteristic (n+1/2) quantum Hall plateaux 
that are found only in single-layer graphene. 
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Supplementary Figure 4: Symmetry in conductance under field reversal B *-* —B, as required 
by reciprocity relations. 
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Supplementary Figure 5: (a) Extracted fluctuations used in Figs. 2 and 3C. Electronic version 
may be zoomed to show full details, (b) An offset and gradient in B± for each value of Vg were 
subtracted from the raw conductance data in order to obtain (a), (c) Variance of a, calculated 
within a sliding window of 0.75V. 
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Supplementary Figure 6: (a) Offsets VQ ffset extracted from Fig. 3c autocorrelations, along 
with autocorrelation side-peak heights relative to the variance. Expected side-peak height is 50% 
for ideal spin-splitting, (b) Fermi speeds vp calculated from VQ ffset . Determination of vf depends 
on the charge neutrality point Vq. Several possible values of Vq are indicated with different colors. 
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Supplementary Figure 7: (a) Duplicated Fig. lg from the paper. Spin-splitting offset increases 
with magnetic field, (b) Idealized extrapolation upwards from the -B||=0 data, showing expected 
appearance of the data if fluctuations split into two identical spin-split copies without otherwise 
changing. 
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Supplementary Figure 8: (a) Autocorrelations for a sliding window (similar to Fig. 3c), with 
a window width of 0.2 volts, and focusing on the region around the Dirac point. A side-peak near 
^offset ^ Q.04V is visible throughout the region near charge neutrality (Vg=0.5-1.5V), correspond- 
ing to a density of states 4 x 10 12 cm _2 /eV. (b) Spin splitting of features with B» near charge 
neutrality. The dashed lines correspond to the density of states extracted from (a). The V's in 
the raw data appear to follow this slope, confirming the constant density of states near charge 
neutrality. 
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Supplementary Figure 9: Negative magnetoconductance due to the in-plane field Bu , expressed 
as a fraction of the zero-field conductance. The curves here were averaged in Vg over 3V-wide 
windows centered at (from bottom to top) -13V,-10V,-7V,-4V,-1V,2V,5V. (5_L=80mT) 
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